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5 ' Abstract 

^ ■ Let S" be a set of n points in M'*. A Steiner convex partition is a tiling of conv(S') with empty 

1^ ■ convex bodies. For every integer d, we show that S admits a Steiner convex partition with at 

^0 . most \{n — l)/d] tiles. This bound is the best possible for points in general position in the 

' plane, and it is best possible apart from constant factors in every fixed dimension d > 3. We 

also give the first constant-factor approximation algorithm for computing a minimum Steiner 
convex partition of a planar point set in general position. 

Establishing a tight lower bound for the maximum volume of a tile in a Steiner convex 
partition of any n points in the unit cube is equivalent to a famous problem of Danzer and 
c/2 . Rogers. It is conjectured that the volume of the largest tile is w(l/n). Here we give a (1 — e)- 

approximation algorithm for computing the maximum volume of an empty convex body amidst 
n given points in the d-dimensional unit box [0, l]''. 

■ 

^ ■ Keywords: Steiner convex partition, Horton set, epsilon-net, lattice polytope, approximation 

' algorithm. 
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1 Introduction 



Let S be a set of n > d + 1 points in R*^, d > 2. A convex body C is empty if its interior is 
disjoint from S. A convex partition of 5" is a partition of the convex hull conv(5') into empty 
convex bodies (called tiles) such that the vertices of the tiles are in S. In a Steiner convex partition 
of S the vertices of the tiles are arbitrary: they can be points in S or Steiner points. For instance, 
^ I any triangulation of 5 is a convex partitions of S, where the convex bodies are simplices, and so 

■ conv(S) can be always partitioned into 0(rel-'^/^-l) empty convex tiles |13]. 

In this paper, we study the minimum number of tiles that a Steiner convex partition of every n 
points in admits, and the maximum volume of a single tile for a given point set. The research 
is motivated by a longstanding open problem by Danzer and Rogers [21 [6l [9l [T8l [33] : What is the 
maximum volume of an empty convex body C C [0, 1]"' that can be found amidst any set S C [0, 1]*^ 
of n points in a unit cube? The current best bounds are 0(l/n) and 0(logn/n), respectively (for 
a fixed d). The lower bound comes from decomposing the unit cube by n parallel hyperplanes, 
each containing at least one point, into at most n + 1 empty convex bodies. The upper bound is 



* A preliminary version of this paper appeared in the Proceedings of the 13th Scandinavian Symposium and Work- 
shops on Algorithm Theory, Helsinki, Finland, July, 2012. 

^Department of Computer Science, University of Wisconsin-Milwaukee, WI 53201-0784, USA. 
Email: duiiiitres@uwin.edu. Supported in part by NSF grant DMS-1001667. 

■'"Department of Computer Science, University of Illinois at Urbana-Champaign, Urbana, IL 61801-2302, USA. 
Email: sariel(3cs .uiuc. edu. Work on this paper was partially supported by a NSF AF award CCF-0915984. 

^Department of Mathematics and Statistics, University of Calgary, Canada; and Department of Computer Science, 
Tufts University, Medford, MA. Email: cdtothOucalgary . ca. Supported in part by NSERC grant RGPIN 35586. 



1 



tight apart from constant factors for n uniformly distributed random points in the unit cube. It is 
suspected that the largest volume is larger than i}{l/n), i.e., it is u}{l)/n in any dimension d>2. 

Minimum number of tiles in a convex partition. A minimum convex partition of 5 is a 
convex partition of S with a minimum number of tiles. Denote this number by fdiS). Further 
define (by slightly abusing notation) 

/rf(n) = max{/rf(5) : 5 C MM^I = n}. 

Similarly define a minimum Steiner convex partition of S as one with a minimum number of tiles 
and let gd{S) denote this number. We also define 

gain) = max{gdiS) : S C \S\ = n}. 

There has been substantial work on estimating f2{n), and computing f2{S) for a given set S 
in the plane. It has been shown successively that f2{n) < 10"--18 j^y Neumann-Lara et al. |32j . 
/2(n) < by Knauer and Spihner i28j, and /2(n) < ^ forn > 6 by Sakai and Urrutia [M]. 

From the other direction, Garci'a-Lopez and Nicolas [19] proved that f2{n) > l^i^z^^ n > 4, 
thereby improving an earlier lower bound f2in) > n + 2 by Aichholzer and Krasser pj. Knauer 
and Spillner [28] have also obtained a jy-factor approximation algorithm for computing a minimum 
convex partition for a given set 5 C M-^, no three of which are collinear. There are also a few exact 
algorithms, including three fixed-parameter algorithms [17 ^ 120 1 [36]. 

The state of affairs is much different in regard to Steiner convex partitions. As pointed out 
in |15] . no corresponding results are known for the variant with Steiner points. Here we take the 
first steps in this direction, and obtain the following results. 

Theorem 1. For n > d + 1, we have gdin) < 1"^^^] • F'or d = 2, this hound is the best possible, 
that is, g2{n) = \{n — l)/2]; and for every fixed d>2, we have gdin) > r2(n). 

We say that a set of points in is in general position if every /c-dimensional affine subspace 
contains at most k + 1 points for < A; < d. We show that in the plane every Steiner convex 
partition for n points in general position, i of which lie in the interior of the convex hull, has at 
least il.{i) tiles. This leads to a simple constant-factor approximation algorithm. 

Theorem 2. Given a set S of n points in general position in the plane, a ratio 3 approximation 
of a minimum Steiner convex partition of S can be computed in O(nlogn) time. 

The average volume of a tile in a Steiner convex partition of n points in the unit cube [0, l]'^ is 
an obvious lower bound for the maximum possible volume of a tile, and for the maximum volume 
of any empty convex body C C [0, 1]"^. The lower bound gdin) > Q{n) in Theorem [1] shows that 
the average volume of a tile is 0(l/n) in some instances, where the constant of proportionality 
depends only on the dimension. This implies that a simple "averaging" argument is not a viable 
avenue for finding a solution to the problem of Danzer and Rogers. 

Maximum empty polytope among n points in a unit cube. In the second part of the 
paper, we consider the following problem: Given a set of n points in rectangular box B in W^, find 
a maximum- volume empty convex body C C B. Since the ratio between volumes is invariant under 
affine transformations, we may assume without loss of generality that B = [0, 1]"^. We therefore have 
the problem of computing a maximum volume empty convex body C C [0, 1]'^ for a set of n points 
m [0,1]'^. It can be argued that the maximum volume empty convex body is a polytope, however. 
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the (number and) location of its vertices is unknown and this represents the main difficulty. For 
d = 2 there is a polynomial-time exact algorithm (see Section Ej) while for d > 3 we are not aware 
of any exact algorithm. Thus the problem of finding faster approximations naturally suggests itself. 

There exist exact algorithms for some related problems. Eppstein et al. [16] find the maximum 
area empty convex /c-gon with vertices among n points in 0{k'n?) time, if it exists. As a byproduct, 
a maximum area empty convex polygon with vertices among n given points can be computed 
exactly in O(n^) time with their dynamic programming algorithm. 

By John's ellipsoid theorem [30], the maximum volume empty ellipsoid in [0,1]°' gives a l/d!^- 
approximation. Here we present a (1 — e)-approximation for the maximum volume empty convex 
body Copt by first guessing a good approximation of the bounding hyperrectangle of Copt of mini- 
mum volume, and then finding a sufficiently close approximation of Copt inside it. We obtain the 
following two approximation algorithms. The planar algorithm runs in near-quadratic time in n, 
however, the running time degrades with the dimension. 

Theorem 3. Given a set S of n points in [0, 1]^ and parameter e > 0, one can compute an empty 
convex body C C [0, 1]^, such that vol(C) > (1 — e)volopt('S'). The running time of the algorithm is 



Theorem 4. Given a set S of n points in [0, 1]*^, d > 3, and a parameter e > 0, one can compute 
an empty convex body C C [0, 1]*^, such that vol(C) > (1 — e)volopt(5'). The running time of the 
algorithm is 



As far as the problem of Danzer and Rogers is concerned, one need not consider convex sets — it 
suffices to consider triangles — and for triangles the problems considered are much simpler. Specif- 
ically, the largest area triangle amidst n points in the unit square U yields a constant-factor ap- 
proximation of the largest area convex body (polygon) amidst the same n points. That is, the 
asymptotic dependencies on n of the areas of the largest empty triangle and convex body are the 
same. Even in higher dimensions, if d is fixed, the asymptotic dependencies on n of the volumes 
of the largest empty simplex and convex body are the same. For d = 2 there is a polynomial-time 
exact algorithm for computing the largest empty triangle contained in U (see Section [5|) while for 
d > 3 we are not aware of any exact algorithm. 

Related work. Decomposing polygonal domains into convex sub-polygons has been also studied 
extensively. We refer to the article by Keil [25] for a survey of results up to the year 2000. For 
instance, when the polygon may contain holes, obtaining a minimum convex partition is NP- 
hard, regardless of whether Steiner points are allowed. For polygons without holes, Chazelle and 
Dobkin [11] obtained an 0(n + r^) time algorithm for the problem of decomposing a polygon with 
n vertices, r of which are reflex, into convex parts, with Steiner points permitted. Keil [25j notes 
that although there are an infinite number of possible locations for the Steiner points, a dynamic 
programming approach is amenable to obtain an exact (optimal) solution; see also |26| I35j. 

Fevens et al. [TT] designed a polynomial time algorithm for computing a minimum convex 
partition for a given set of n points in the plane if the points are arranged on a constant number 
of convex layers. The problem of minimizing the total Euclidean length of the edges of a convex 
partition has been also considered. Grantson and Levcopoulos [19j, and Spillner [36j proved that 
the shortest convex partition and Steiner convex partition problems are fixed parameter tractable, 
where the parameter is the number of points of P lying in the interior of conv(P). Dumitrescu 
and Toth [T5] proved that every set of n points in admits a Steiner convex partition which is at 
most 0(logn/loglogn) times longer than the minimum spanning tree, and this bound cannot be 
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improved. Without Steiner points, the best upper bound for the ratio of the minimum length of a 
convex partition and the length of a minimum spanning tree (MST) is 0{n) j27j . 

A largest area convex polygon contained in a given (non-convex) polygon with n vertices can be 
found by the algorithm of Chang and Yap [10] in 0{n'^) time. The problem is known as the potato- 
peeling problem. On the other hand, a largest area triangle contained in a simple polygon with 
n vertices, can be found by the algorithm of Melissaratos and Souvaine [STj in 0{n'^) time. Hall- 
Holt et al. J2T\ compute a constant approximation in time 0(n log n). The same authors show how 
to compute a (1 — e)-approximation of the largest fat triangle inside a simple polygon (if it exists) 
in time 0{n). Given a triangulated polygon (with possible holes) with n vertices, Aronov et al. [1] 
compute the largest area convex polygon respecting the triangulation edges in 0{v?) time. 

For finding a maximum volume empty axis-parallel box amidst n points in [0, l]'^. Backer and 
Keil [5J reported an algorithm with worst-case running time of 0{n'^\og'^~'^ n). An empty axis- 
aligned box whose volume is at least (1— e) of the maximum can be computed in O (^(^)'^ \og'^ 
time by the algorithm of Dumitrescu and Jiang [H] . 

Lawrence and Morris [29] studied the minimum integer kd{n) such that the complement M'^\S' of 
any n-element set S C W^, not all in a hyperplane, can be covered by kd{n) convex sets. They prove 
kd{n) > 0(log n/dlog log n). It is known that covering the complement of n uniformly distributed 
points in [0, 1]"^ requires at least ft{n/dlogn) convex sets, which follows from the upper bound in 
the problem of Danzer and Rogers. 

2 Combinatorial bounds 

In this section we prove Theorem[TJ We start with the upper bound. The following simple algorithm 
returns a Steiner convex partition with at most \{n — l)/d'\ tiles for any n points in M'^. 

Algorithm Al: 

Step 1. Compute the convex hull R ^ conv(S') of S. Let ^4 C be the set of hull vertices, and 
let B = S \A denote the remaining points. 

Step 2. Compute conv(i?), and let H be the supporting hyperplane of an arbitrary {d — 1)- 
dimensional face of conv(i?). Denote by the halfspace that contains B, and H~ = M.'^\H^. 
The hyperplane H contains d points of B, and it decomposes R into two convex bodies: RriH~ 
is empty and R ^ ROH^ contains all points in B\H. Update B ^ B\H and R ^ RCiH^. 

Step 3. Repeat Step 2 with the new values of R and B until B is the empty set. (If \B\ < d, 
then any supporting hyperplane of B completes the partition.) 




Figure 1: Steiner convex partitions with Steiner points drawn as hollow circles. Left: A Steiner convex 
partition of a set of 13 points. Middle: A Steiner partition of a set of 12 points into three tiles. Right: A 
Steiner partition of the same set of 12 points into 4 tiles, generated by Algorithm Al (the labels reflect the 
order of execution). 



4 



It is obvious that the algorithm generates a Steiner convex partition of S. An illustration 
of Algorithm Al on a small planar example appears in Figure [U (right). Let h and i denote 
the number of hull and interior points of S, respectively, so that n = h + i. Each hyperplane 
used by the algorithm removes d interior points of S (with the possible exception of the last 
round if i is not a multiple of d). Hence the number of convex tiles is 1 + [i/d], and we have 
1 + \i/d~\ = \{i + d)/d~\ < \{n - l)/d~\, as required. 

Lower bound in the plane. A matching lower bound in the plane is given by the following 
construction. For n > 3, let S* = A U where ^ is a set of 3 non-collinear points in the plane, 
and i3 is a set of n — 3 points that form a regular (n — 3)-gon in the interior of conv(^), so that 
conv(5) = conv(A) is a triangle. If n = 3, then conv(5') is an empty triangle, and g2{S) = 1. If 
4 < Ti < 5, is not in convex position, and so g2{S) > 2. Suppose now that n > 6. 

Consider an arbitrary convex partition of S. Let o be a point in the interior of conv(i3) such 
that the lines os, s G S, do not contain any edges of the tiles. Refer to Figure [2] (left). For each 
point s G B, choose a reference point r(s) G M? on the ray ot in conv(j4) \conv(i?) sufficiently close 
to point s, and lying in the interior of a tile. Note that the convex tile containing o cannot contain 
any reference points. We claim that any tile contains at most 2 reference points. This immediately 
implies 52 (5") >l+\{n- 3)/2] = \{n - l)/2]. 

Suppose, to the contrary, that a tile r contains 3 reference points ri,r2,r3, corresponding to 
the points si, 52,83. Refer to Figure [2j Note that o cannot be in the interior of r, otherwise r 




Figure 2: Left: Lower bound construction in M^. Right: Points in general position on a saddle surface in R^. 

would contain all points si,S2,ss in its interior. Hence conv{o, si, S2, S3} is a quadrilateral, and 
convjo, ri, r2, rs} is also a quadrilateral, since the reference points are sufficiently close to the 
corresponding points in B. We may assume w.l.o.g. that vertices of conv{o, si, S2, S3} are o, si, S2, 
S3 in counterclockwise order. Then S2 lies in the interior of conv{o, ri, r2, r3}. We conclude that 
every tile r contains at most 2 reference points, as required. 

Lower bounds for d > 3. A similar construction works in for any d > 2, but the lower bound 
no longer matches the upper bound gdin) < \{n — l)/d] for d > 3. 

Recall that a Norton set |24] is a set S oi n points in the plane such that the convex hull of any 
7 points is non-empty. Valtr [37] generalized Horton sets to M*^. For every d € N, there exists a 
minimal integer h{d) with the property that for every n G N there is a set of n points in general 
position in such that the convex hull of any h{d) + 1 points in S is non-empty. It is known that 
h{2) = 6, and Valtr proved that h{3) < 22, and in general that h{d) < 2'^~'^{N{d - 1) + 1), where 
N{k) is the product of the first k primes. 
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We construct a set 5 of n > d + 1 points in R*^ as follows. Let S = Au B, where ^ is a set of 
d+1 points in general position in W^, and 5 is a generalized Horton set of n — (d+ 1) points in the 
interior of conv(^), such that the interior of any h{d) + 1 points from B contains some point in B. 

Consider an arbitrary Steiner convex partition of S. Every point 6 G i? is in the interior of 
conv(5), and so it lies on the boundary of at least 2 convex tiles. For each b & B, place two reference 
points in the interiors of 2 distinct tiles incident to b. Every tile contains at most h{d) reference 
points. Indeed, if a tile contains h{d) + 1 reference points, then it is incident to h{d) + 1 points in 
B, and some point of B lies in the interior of the convex hull of these points, a contradiction. 

We have 2(n — d — 1) reference points, and every tile contains at most h{d) of them. So the 
number of tiles is at least \2{n — d — l)/h{d)'\ . 



3 Approximating the minimum Steiner convex partition in 

In this section we prove Theorem [2] by showing that our simple-minded algorithm Al from Section [2] 
achieves a constant-factor approximation in the plane if the points in S are in general position. 



Approximation ratio. Recall that algorithm Al computes a Steiner convex partition of conv(S') 
into at most 1 + \i/2\ parts, where i stands for the number of interior points of S. 

If i = 0, the algorithm computes an optimal partition, i.e., ALG = OPT = 1. Assume now that 

1 > 1. Consider an optimal Steiner convex partition 11 of 5 with OPT tiles. We construct a planar 
multigraph G = {V, E) as follows. The faces of G are the convex tiles and the exterior of conv(5') 
(the outer face). The vertices V are the points in the plane incident to at least 3 faces (counting 
the outer face as well). Since i > 1, G is non-empty and we have \V\ > 2. Each edge in is a 
Jordan arc on the common boundary of two faces. An edge between two bounded faces is a straight 
line segment, and so it contains at most two interior points of S. An edge between the outer face 
and a bounded face is a convex arc, containing hull points from S. Double edges are possible if 
two vertices of the outer face are connected by a straight line edge and a curve edge along the 
boundary — in this case these two parallel edges bound a convex face. No loops are possible in G. 
Since 11 is a convex partition, G is connected. 

Let V, e, and /, respectively, denote the number of vertices, edges, and bounded (convex) faces 
of G; in particular, / = OPT. By Euler's formula for planar multigraphs, we have v — e + f = 1, 
that is, f = e — V + 1. By construction, each vertex of G is incident to at least 3 edges, and every 
edge is incident to two vertices. Therefore, 3v < 2e, or v < 2e/3. Consequently, f = e — v + 1 > 
e — 2e/3 + l = e/3 + 1. Since S is in general position, each straight-line edge of G contains at most 

2 interior points from S. Curve edges along the boundary do not contain interior points. Hence 
each edge in E is incident to at most two interior points in S, thus i < 2e. Substituting this into 
the previous inequality on / yields OPT = / > e/3 + 1 > z/6 + 1. Comparing this lower bound 
with the upper bound ALG < [i/2] + 1, we conclude that 

ALG ^ \i/2]+l ^ i + 3 
OPT - i/6 + 1 - i + 6 

and the approximation ratio of 3 follows. 



Tightness of the approximation ratio. We first show that the above ratio 3 is tight for 
Algorithm Al. We construct a planar point set S as follows. Consider a large (say, hexagonal) 
section of a hexagonal lattice. Place Steiner vertices at the lattice points, and place two points 
in S on each lattice edge. Slightly perturb the lattice, and add a few more points in S near the 
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boundary, and a few more Steiner points, so as to obtain a Steiner convex partition of S with no 
three points colhnear. Denote by v, e, and /, the elements of the planar multigraph G as before. 
Since we consider a large lattice section, we have v,e, f — >• oo. We write a ^ b, whenever a/b — t- 1. 
As before, we have f + v = e + 1, and since each non-boundary edge is shared by two convex 
faces, we have e ~ 6//2 = 3/. By construction, i ~ 2e ~ 6/, hence / ~ i/6. Therefore the 
convex partition constructed above has / ~ i/6, while Algorithm Al constructs one with about 
i/2 faces. Letting e — )• oo, then i — t- oo, and the ratio ALG/OPT approaches 3 in the limit: 
ALG/OPT ~ (i/2)/(i/6) = 3. 

Time analysis. It is easy to show that Algorithm Al runs in 0(n log n) time for a set S of 
n points in the plane. We employ the semi-dynamic (delete only) convex hull data structure of 
Hershberger and Suri [23]. This data structure supports point deletion in O(logn) time, and uses 
0(n) space and O(nlogn) preprocessing time. We maintain the boundary of a convex polygon R 
in a binary search tree, a set B G S oi points lying in the interior of R, and the convex hull conv(i?) 
with the above semi-dynamic data structure [23]. Initially, R = conv(5'), which can be computed 
in 0(n log n) time; and B C S is the set of interior points. In each round of the algorithm, consider 
the supporting line H of an arbitrary edge e of conv(i?) such that B lies in the halfplane . The 
two intersection points of H with the boundary of R can be computed in O(logn) time. At the end 
of the round, we can update B ^ B\H and conv(S) in 0{k\ogn) time, where k is the number of 
points removed from B] and we can update R ^ RCiH'^ in O(logn) time. Every point is removed 
from B exactly once, and the number of rounds is at most \{n — 3)/2] , so the total update time is 
0(n log n) throughout the algorithm. 

Remark. Interestingly enough, in dimensions 3 and higher. Algorithm Al does not give a 
constant-factor approximation. For every integer n, one can construct a set 5 of n points in 
general position in such that i = n — 4 of them lie in the interior of conv(5), but the minimum 
Steiner convex partition has only 0(-y/n) tiles. In contrast. Algorithm Al computes a Steiner 
partition with i/3 = (n — 4)/3 convex tiles. 

We first construct the convex tiles, and then describe the point set S. Specifically, S consists 
of 4 points of a large tetrahedron, and 3 points in general position on the common boundary of 
certain pairs of adjacent tiles. 

Let k = \y/ {n — 4)/3] . Place {k + 1)^ Steiner points {a,h,a^ — 6^) on the saddle surface 
z = — for pairs of integers (a, 6) G Z^, — [/c/2j < a,b < \k/2\. The four points {(x, y, — y^) : 
X G {a, a + 1}, y € {b,b+ 1}} form a parallelogram for every (a, h) G 1? , — \k/2\ < a,b < \k/2] — 1. 
Refer to Figure [2] (right). These parallelograms form a terrain over the region {(x,y) : —[k/2\ < 
x,y < \k/2\}. Note that no two parallelograms are coplanar. Subdivide the space below this 
terrain by vertical planes x = a, —[k/2\ < a < [fc/2]. Similarly, subdivide the space above this 
terrain by planes y = b, —[k/2\ < b < \k/2'\. We obtain 2k interior-disjoint convex regions, k 
above and k below the terrain, such that the common boundary of a region above and a region 
below is a parallelogram of the terrain. The points in that do not lie above or below the terrain 
can be covered by 4 convex wedges. 

Enclose the terrain in a sufficiently large tetrahedron T. Clip the 2k convex regions and the 4 
wedges into the interior of T. These 2A; + 4 convex bodies tile T. Choose 3 noncollinear points of 
S in each of the k'^ parallelograms, such that no 4 points are coplanar and no 2 are collinear with 
vertices of T. Let the point set S be the set of 4 vertices of the large tetrahedron T and the 3/c^ 
points selected from the parallelograms. 
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4 Approximating the maximum empty convex body 



Let S be a set of points in the unit cube [0, l]'^ C M'^. Our task is to approximate the largest convex 
body C C [0, l]'' that contains no points of S in its interior. Let Copt = Copt{S) denote this body, 
and let volopt(5') denote its volume. 



4.1 Approximation by the discrete hull 

In the following, assume that m > is some integer, and consider the grid point set 

Q{m) = . . .,id)/m ii, . . . , G {0, 1, . . . , m} | . 

Let S C [0, l]'^ be a point set, and let Copt be the corresponding largest empty convex body in 
[0, 1]*^. Given a grid Gim), we call conv(Copt n Gim)) the discrete hull of Copt [22]. We need the 
following easy lemma. 

Lemma 1. Let C C [0, l]'^ he a convex body and D = conv(C n ^(m)). Then we have vol(C) — 
vol(L') = 0{l/m), where the constant of proportionality depends only on d. 

Proof. Consider a point p £ C \ D. Consider the set of 2d points X = {p± 2{d/m)ei}, where e, 
is the unit vector having one in the ith coordinate, and everywhere else. If one of the points of 
X is outside C, then the distance from p to the boundary of C is at most 2d/m. Otherwise, the 
cube p+ [-2,2]"^/ m is contained in the "diamond" conv(X), which is in turn contained in C. But 
then, the grid points of the grid cell of G{rn) containing p are in C, and p can not be outside D. 
We reached a contradiction. 

It follows that all the points in the corridor C \ D are at distance at most 2{d/m) from the 
boundary of C. The volume of the boundary of C is bounded by the volume of the boundary of 
the unit cube, namely 2d. As such, the volume of this corridor is vol(9C) 0{d/m) < {2d){2d/m) = 
0{d^ /m). For a fixed d, this is 0(l/m), as claimed. □ 

Lemma [1] implies that if volopt('S') > />, in order to obtain a (1 — e)-approximation, we can 
concentrate our search on convex polytopes that have their vertices at grid points in G{m), where 
m = 0{l/{ep)). If p is a constant, then the maximum volume empty lattice polytope in ^(m) with 
m = 0{l/e) is an (1 — e)-approximation for Copt- However, for arbitrary voloptC-S) > 0(l/n), a 
much finer grid would be necessary to achieve this approximation. 



4.2 An initial brute force approach 

In this section we present approximation algorithms relying on Lemma [1] alone, approximating the 
maximum volume empty polytope by a lattice polytope in a sufficiently fine lattice (grid). We shall 
refine our technique in Subsections 14.31 and 14.41 

For the plane, we take advantage of the existence of an efficient solution for a related search 
problem in the spirit of those considered by Eppstein et al. |16| . 

Lemma 2. Given a set S of n points and a set Q of m points in the plane, one can compute a 
convex polygon of the largest area with vertices in S that does not contain any point of Q in its 
interior in 0(n'^m + n^) time. The algorithm has the same running time ifQ is a set ofm forbidden 
rectangles. 
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Proof. The algorithm works by dynamic programming. First, we compute for all triangles with 
vertices from S whether they contain a forbidden point inside them; trivially this can be done in 
0{rfim) time. We then build a directed graph G on the allowable triangles, connecting two triangles 
A and A' if they share their left endpoint, are interior disjoint, share an edge, and their union forms 
a convex quadrilateral. We orient the edge from the triangle that is most counterclockwise (around 
the common vertex) to the other triangle. All edges are oriented "upwards", so G is a directed 
acyclic graph (DAG). Observe that G has 0{rfi) vertices (allowable triangles) and the maximum 
out-degree in G is bounded from above by n. 

The weight of a vertex corresponding to a triangle is equal to its area. Clearly, a convex polygon 
corresponds to a path in G, namely the triangulation of the polygon from its leftmost vertex, and 
its weight is the area of the polygon. Finding the maximum weight path can be done in linear time 
in the size of the DAG; see e.g., [IH Section 4.7]. G has 0{v?) vertices and O(n^) edges, and as 
such the overall running time is 0[n^m + n'^). □ 

The combination of Lemmas [T] and [2] readily yields an approximation algorithm for the plane, 
whose running time depends on vol ( Copt )• 

Lemma 3. Given a set S C [0, 1]^ ofn points, such that volopt('S') > p, and a parameter e > 0, one 
can compute an empty convex body C C [0,1]^ such that vol(C) > (1 — e)volopt(>S'). The running 
time of the algorithm is 0(n + l/(ep)®) . 

Proof. Consider the grid G{m) with m = 0(l/(ep)). By Lemma[I]we can restrict our search to a 
grid polygon. Going a step further, we mark all the grid cells containing points of S as forbidden. 
Arguing as in Lemma [H one can show that the area of the largest convex grid polygon avoiding 
the forbidden cells is at least volopt(5') — c/m, where c is a constant. 

We now restrict our attention to the task of finding a largest polygon. We have a set Q of 
0{m?) grid points that might be used as vertices of the grid polygon, and a set of 0{m?) grid cells 
that can not intersect the interior of the computed polygon. By Lemma O a largest empty polygon 
can be found in 0{m^) time. Setting m = 0(1/ (ep)), we get an algorithm with overall running 
time 0{n + l/{ep)^). □ 

In dimensions d > 3, we are not aware of any analogue of the dynamic programming algorithm 
in Lemma [2J Instead, we fall back on a brute force approach that enumerates all feasible subsets 
of a sufficiently fine grid. 

Lemma 4. Given a set S Q [0, 1]^ ofn points, such that volopt(<S') > p, and a parameter e > 0, one 
can compute an empty convex body C C [0, 1]*^, such that vol(C) > (1 — e)volopt('S'). The running 
time of the algorithm is 0{n) + exp(0(m'^('^^"'^^/*^'^+^^ logm)), where m = 0{l/{ep)) and d is fixed. 

Proof. Consider the grid Q{m) with m = 0{1/ {ep)). Let X be the set of vertices of all grid cells 
of Q{m) that contain some point from S (i.e., 2*^ vertices per cell). Note that \X\ = 0{m'^). 
Andrews ^ proved that a convex lattice polytope of volume V has 0(l/('^~^)/('^+^)) vertices. Hence 
a convex lattice polytope in Q{m) has 0{rn'^^'^^'^'>/'^'^^^^) vertices. By the well-known inequality 
(fc) < {^f, the number of subsets of size 0(m'^{'^-i)/(rf+i)) from g{ni) is 

(o(™<"l)/Wi),) ^ (™""'"") ' <exp(o(™<<'-WWi),og,„)) , 

For each such candidate subset G of size 0{rn'^^'^~^^/^'^^'^'^), test whether conv(G) is empty of points 
from X. For each point in X, the containment test reduces to a linear program that can be solved 
in time polynomial in m. Returning the subset with the largest hull volume found yields the desired 
approximation. The runtime of the algorithm is exp 
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Barany and Vershik [8] proved that there are exp (^0{m'^^'^~^^/^'^~^^^)^ convex lattice polytopes 
in Q{m). If the polytopes can also be enumerated in this time (as in the planar case [7]), then the 
runtime in Lemma H] reduces accordingly. 

4.3 A better approximation in the plane 

If Copt is long and skinny (e.g., p is close to 1/n), then the uniform grid G{m) we used in LemmasO 
and m is unsuitable for finding a (1 — e)-approximation efficiently. Instead, we employ a rotated 
and stretched grid (an afhne copy of G{m)) that has similar orientation and aspect ratio as Copt- 
This overcomes one of the main difficulties in obtaining a good approximation. Since we do not 
know the shape and orientation of Copt , we guess these parameters via the minimum area rectangle 
containing Copt- 

Lemma 5. Given a set S C [0, 1]^ of n points with volopt('S') > p and a parameter e > 0, one can 
compute an empty convex body C C [0, 1]^ such that vol(C) > (1 — e)volopt('S'). The running time 
of the algorithm is 0(p~^e~^ nlog^ n) . 

Proof. The intuitive idea is now to guess a rectangle R that contains Copt such that its area, 
vol(Copt), is at least a constant fraction of the area of -R, namely vol(-R); and then apply Lemma[3] 
to the rectangle R (as the unit square) to get the desired approximation. 

Let Bq be the minimum area rectangle (or arbitrary orientation) that contains Copt; see Fig- 
ure [3] (left). We guess an approximate copy of Bq. In particular, we guess the lengths of the two 
sides of Bq (up to a factor of 2) and the orientation of Bq (up to an angle of 0(l/n)), and then try 
to position a scaled copy of the guessed rectangle so that that it fully contains Copt- 





Figure 3: Left: Copt, a niinimuni area rectangle B^, Copt ^ B^, and a minimum area rectangle Bi, Bq (- Bi, 
with canonical side lengths and the same orientation as Bq. Right: Rectangle Bi, a rotated copy B2 with 
the closest canonical orientation, and a minimum area scaled copy kB2 such that Bi C fci?2. 



Assume for convenience that n > 10. Denote by a and b the lengths of the two sides of rectangle 
Bq, where a < b. It is clear that b < \/2, the diameter of the unit square. We now show that 



n-1+4- 



> 



-i-T > = ^, for n > 



10. 



a > 1/n. From Theorem [H we have vol(Copt) > l/[ 

From the other direction, vol(Copt) < \/2 width(Copt) < \/2a. Comparing these bounds yields 
a > 1/n. It follows that if we want to guess the side lengths of Bq up to a factor of two, we need 



to consider only 0(log n) possibilities. Indeed, we consider the canonical lengths 



2Yn, for 



10 



i = 0, . . . , [log2 n\ . Let Bi be a minimum area rectangle with canonical side lengths and the same 
orientation as Bq, so that Bq <^ Bi. 

The orientation of a rectangle is given by the angle between one side and the x-axis. We 
approximate the orientation of Bq by the n canonical orientations a = iir/ (5n), for i = 0,1, . . . , 5n— 
1. Let B2 be a congruent copy of Bi rotated clockwise to the nearest canonical orientation about 
the center of Bi. We show that Bi C 2B2, i.e., a scaled copy of B2 contains Bi. Let k > 1 
be the minimum scale factor such that Bi C kB2- Refer to Figure [3] (right). Denote by the 
common center of Bi and B2, let x be a vertex of Bi on the boundary of kB2, and let y be the 
corresponding vertex of kB2. Clearly, sm(Zxoy) < 7r/(5n) since we rotate by at most 7r/(5n). We 
also have sin(Zoyx) > (l/n)/\/2 > 7r/(5n), since sin(Zoyx) is the ratio between a side of kB2 and 
its diameter (and this ratio is the same for Bq). The law of sines yields \ox\ > \xy\; and we have 
I ox I + I XT/ 1 > \oy\ by the triangle inequality. If follows that \oy\ < 2\ox\, and so A; < 2 suffices. 

So far we have guessed the canonical side lengths and orientation of B2, however, we do not 
know its location in the plane. If a translated copy B2 + v of B2 intersects Copt, then 3B2 + v 
contains it, since Copt ^ -Bq C i?i C 2B2. Consider an arbitrary tiling of the plane with translates 
of B2. By a packing argument, only 0{l/p) translates intersect the unit square [0, 1]^. One of these 
translates, say B2 + v, intersects Copt, and hence the rectangle R = 3B2 + v contains Copt. 

We can apply Lemma[3]to the rectangle R (as the unit square) to get the desired approximation. 
Specifically, let T : — )• be an affine transformation that maps R into the unit square [0, 1]^, 
and apply Lemma [3] for the point set T{S (1 R) and T(R n [0, 1]^). The grid Q{m) clipped in 
T{Rn [0, 1]^) corresponds to a stretched and rotated grid in R; each grid cell of Q{m) is stretched 
to a rectangle with the same aspect ratio as R. The convex polygon Copt occupies a constant 
fraction of the area of R, and so the resulting running time is 0(l/e^). Note that the algorithm 
of Lemma [3] partitions R into a grid with 0(l/e^) cells. The approximation algorithm only cares 
about which cells are empty and which are not. 

Since we have to repeat this for all possible side lengths, orientations, and positions of R, the 
overall running time is 




Since p = r2(l/n), the above expression is bounded by 0(e ^n? log^ nj, as claimed. □ 

Remark. If /? = ^^(1) the running time of this planar algorithm is near linear in n. 

By doing an exponentially decreasing search for p, the running time increases only by a constant 
factor (this is a geometrically decreasing series, hence the term with the last value of p dominates 
the whole running time). We summarize our result for the plane in the following. 

Theorem [S], Given a set S of n points in [0, 1]^ and a parameter e > 0, one can compute an empty 
convex body C C [0, 1]^, such that vol(C) > (1 — e)volopt(S'). The running time of the algorithm is 
0[e~^v? log^ n) . 

4.4 A better approximation in higher dimensions 

Given a set 5 C [0, l]'^ of n points and a parameter e > 0, we compute an empty convex body 
C C [0, 1]^ such that vol(C) > (1 — e)volopt('S'). Similarly to the algorithm in Subsection 14.31 we 
guess a hyperrectangle R that contains Copt such that vol(Copt) is at least a constant fraction of 
vol(i?); and then apply Lemma H] to R (as the hypercube) to obtain the desired approximation. 
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Consider a hyperrectangle Bq of minimum volume (and arbitrary orientation) that contains 
Copt- The d edges incident to a vertex of a hyperrectangle B are pairwise orthogonal. We call these 
d directions the axes of B; and the orientation of B is the set of its axes. 

We next enumerate all possible discretized hyperrectangles of volume ^}{l/n), guessing the 
lengths of their axes, their orientations, and their locations as follows: 

Guess the length of every axis up to a factor of 2. Since the minimum length of an axis in 
our case is J7(l/n) and the maximum is Vd, the number of possible lengths to be considered is 
0(log°'n). Let Bi be a hyperrectangle of minimum volume with canonical side lengths and the 
same orientation as Bq such that Bq Bi. 

We can discretize the orientation of a hyperrectangle as follows. We spread a dense set of points 
on the sphere of directions, with angular distance 0(l/n) between any point on the sphere and its 
closest point in the chosen set. 0(n^~^) points suffice for this purpose. We try each point as the 
direction of the first axis of the hyperrectangle, and then generate the directions of the remaining 
axes analogously in the orthogonal hyperplane for the chosen direction. Overall, this generates 
0(nS^=i = 0(n'^('^-i)/2) possibilities. 

Successively replace each axis of Bi by an approximate axis that makes an angle at most 
a = l/(cn) with its corresponding axis, where c = c{d) is a constant depending on d. Let B2 be a 
congruent copy of Bi obtained in this way. If c = c{d) is sufficiently small, then Bi C 2B2- 

Consider a tiling of with translates of B2. Note that only 0(l/vol(Copt)) = 0{n) trans- 
lates intersect the unit cube [0,1]'^. One of these translates B2 + v intersects Copt, and then the 
hyperrectangle R = 3B2 + v contains Copt. Since Copt (5') takes a constant fraction of the vol- 
ume of R, we can deploy Lemma H] in this case, and get the desired (1 — e)-approximation in 
exp(0(e~'^^'^~^)/('^"'"^) loge~^)) time. Putting everything together, we obtain the following. 

Theorem |4] Given a set S of n points in [0,1]°', d > 3, and a parameter e > 0, one can compute 
an empty convex body C C [0, l]'^, such that vol(C) > (1 — e)volopt(>S'). The running time of the 
algorithm is 

exp(o(e-^('^-i)/(^+i) loge-i)) ni+^('^-i)/2 j^gd^^ 

Remark. Consider a set 5 of n points in W^. The approximation algorithm we have presented can 
be modified to approximate the largest empty tile, i.e., the largest empty convex body contained in 
conv(S), rather than [0, l]'^. The running time is slightly worse, since we need to take the boundary 
of conv(5) into account. We omit the details. 

5 Conclusions 

In this section we briefly outline two exact algorithms for finding the largest area empty convex 
polygon and the largest area empty triangle amidst n points in the unit square. At the end we list 
a few open problems. 

Largest area convex polygon. Let S C U = [0,1]^, where |S| = n. Let T be the set of four 
vertices of U. Observe that the boundary of an optimal convex body, Copt, contains at least two 
points from S U T. By convexity, the midpoint of one of these O(n^) segments lies in Copt- For 
each such midpoint m, create a weakly simple polygon Pm by connecting each point p G S to the 
boundary of the square along the ray mp. The polygon has 0{n) vertices and is empty of 
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points from S in its interior. Then apply the algorithm of Chang and Yap [TO] for the potato- 
peeling problem (mentioned in Section [T]) in these O(n^) weakly simple polygons. The algorithm 
computes a largest area empty convex polygon contained in a given (non-convex) polygon with n 
vertices in O(n^) time. Finally, return the largest convex polygon obtained this way. The overall 
running time is 0{n^). 

The running time can be reduced to 0(n^ log n) as follows. Instead of considering the 0{n?) 
midpoints, compute a set P of O(nlogn) points so that every convex set of area at least 2/(n + 4) 
contains at least one of these points. In particular, Copt contains a point from P. The set P 
can be computed by starting with a 0{n) x 0(n) grid, and then computing an e-net for it, where 
e = 0(l/n), using discrepancy |30) . The running time of this deterministic procedure is roughly 
O(n^), and the running time of the overall algorithm improves to 0(n'' • nlogn) = O(n^logn). 

Largest area empty triangle. The same reduction can be used for finding largest area empty 
triangle contained in U, resulting in 0{'n?) weakly simple polygons Pm- Then the algorithm of 
Melissaratos and Souvaine [31j for finding a largest area triangle contained in a polygon is applied 
to each of these O(n^) polygons. The algorithm computes such a triangle in O(n^) time, given a 
polygon with n vertices. Finally, return the largest convex polygon obtained this way. The overall 
running time is O(n^). Via the e-net approach (from the previous paragraph) the running time of 
the algorithm improves to 0{n'^ • nlogn) = 0{n^ logn). 

Open questions. Interesting questions remain open regarding the structure of optimal Steiner 
convex partitions and the computational complexity of computing such partitions. Other questions 
relate to the problem of finding the largest empty convex body in the presence of points. We list 
some of them: 

(1) Is there a polynomial-time algorithm for computing a minimum Steiner convex partition of a 
given set of n points in M*^? Is there one for points in the plane? 

(2) Is there a constant-factor approximation algorithm for the minimum Steiner convex partition 
of an arbitrary point set in M'^ (without the general position restriction)? Is there one for 
points in the plane? 

(3) For d > 2, the running time of our approximation algorithm for the maximum empty polytope 
has a factor of the form rf^'^^\ It seems natural to conjecture that this term can be reduced 
to rP^'^\ Another issue of interest is extending Lemma [2] to higher dimensions for a faster 
overall algorithm. 

(4) Given n points in [0, 1]'^, the problem of finding the largest convex body in [0, 1]"^ that contains 
up to k (outlier) points naturally suggests itself and appears to be also quite challenging. 

Acknowledgement. The authors thank Joe Mitchell for helpful discussions regarding the exact 
algorithms in Section [H in particular for suggesting the reduction of the maximum-area-empty- 
convex-body problem to the potato-peeling problem. 
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